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Abstract. We study theoretically the collective E(g)e Jahn-Teller-Dicke distortion in 
a system of trapped ions. We focus in the limit of infinite range interactions in which an 
ensemble of effective spins interacts with two collective vibrational modes with U(l) 
symmetric couplings. Our model is exactly solvable in the thermodynamical limit 
and it is amenable to be solved by exact numerical diagonalization for a moderate 
number of ions. We show that trapped ions are ideally suited to study the emergence 
of spontaneous symmetry breaking of a continuous symmetry and magnetic structural 
phase transition in a mcsoscopic system. 
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1. Introduction 

Physical systems wliere bosonic modes interact with electronic or pseudospin degrees 
of freedom reveal a rich variety of phenomena in condensed matter and atomic physics. 
A prominent example is given by the Jahn- Teller (JT) models [1,2] which describe the 
interaction of electronic orbital degrees of freedom with vibrational modes either in 
molecules or solids. The JT effect is formulated as a structural instability of molecular 
configurations in electronically degenerate states. In particular, the electron-phonon 
coupling shifts the potential minima of the nuclei, which leads to position reordering 
and molecular distortion. Similar to the molecular systems, the properties of some 
crystals are also strongly affected by the JT coupling, including symmetry breaking and 
structural phase transitions [3]. Furthermore, the strong electron-phonon coupling in 
the cooperative JT models is an important factor in the description of colossal magneto- 
resistance in manganites and high Tc-superconductivity [4,5]. 

Atomic systems such as ultracold atoms and trapped ions allow experimentalists 
to implement JT models in a controllable way that is not possible in solid-state or 
molecular setups. This is a motivation to push the current quantum technology toward 
the realization of Analogical Quantum Simulators (AQS). The latter are controllable 
systems where interactions between particles can be tuned and quantum states can be 
accurately prepared and measured with high efficiency. Recently, physical realizations 
of JT couplings have been discussed in terms of two-level systems coupled to a bimodal 
cavity [6] and Bose-Einstein condensates in the presence of spatially dependent laser 
fields [7]. These systems pave the way for studying a quantum phenomena such as 
ground-state entanglement [8,9] and creation of artificial non-Abelian magnetic fields 
[10]. A quantum chaotic behavior in the energy spectrum of multi-spin lattice JT model 
was discussed in [11]. 

Among the most promising physical systems for implementing AQS are linear ion 
crystals interacting with external lasers or magnetic fields [12, 13]. The main advantages 
of trapped ions are their addressability, long coherence times and high fidelity readout. 
The current available ion trapping technology allows to explore the physics of quantum 
phase transitions in complex spin systems [14-18], interacting bosons [19-23], relativistic 
effects, [24,25] and quantum open systems [26]. 

In this work we propose an implementation of AQS of an infinite range E^e Jahn- 
Teller-Dicke (JTD) model based on a trapped ion crystal. The doublet of electronic 
states is replaced here by two internal states of the ions. The pair of molecular 
vibrational modes is represented by the two degenerate orthogonal center-of-mass (cm) 
modes, respectively, in the two radial directions. We show that the U(l) symmetric 
JTD spin-phonon coupling can be provided by applying a magnetic field with time- 
oscillating gradient, which couples the collective spin ensemble to the two orthogonal 
cm modes [27-29], Fig. 1. Such an oscillating magnetic field has been used to implement 
experimentally a single qubit rotation and multi-qubit quantum gates [30] , which shows 
that our model could be realized with the current state-of-the art ion trap setup. In a 
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Figure 1. Exemplary geometry for the realization of the collective Jahn- Teller-Dicke 
model with time- varying magnetic field gradient produced by a surface-electrode trap. 
A linear ion crystal is oriented along the z axis. Oscillating currents in two wires Wi 
and W2 parallel to the linear ion crystal are used to provide the symmetric spin-phonon 
coupling. 



previous work [21] we have focused on the cooperative JT model in which all vibrational 
modes are coupled to the effective spins, which leads to spin-phonon quasi-condensate. 
Here, we show that under suitable experimental conditions the spin ensemble may 
interact only with a single vibrational mode in each radial direction. Thus, we deal 
with infinite range JT models, which are quasi-exactly solvable, in the sense that the 
ground state can be found in the thermodynamical limit [31] or even studied by exact 
numerical diagonalization with a moderate number of ions. 

The E(8>e JTD model possesses a continuous symmetry associated with rotation in 
the plane orthogonal to the trap axis. There is a critical spin-phonon coupling above 
which the U(l) symmetry is spontaneously broken and the system evolves into one 
particular ground state which does not respect the same symmetry as the Hamiltonian. 
That broken symmetry is associated with a magnetic structural phase transition at zero 
temperature, where ions' equilibrium positions are arbitrarily displaced in the radial 
x-y plane together with the creation of macroscopic spin coherence. We show that the 
radial distortion of the ion crystal increases with the number of ions as a/ZV, while the 
mean phonon number and the spin coherence scale as N, which is an analog to the 
normal-to-super-radiance phase transition in the Dicke model [32,33]. 

The paper is arranged as follows: In Sec. 2 we describe the collective vibrational 
spectrum of the linear ion crystal. We show that an oscillating magnetic field gradient 
could select only one vibrational mode in each radial direction and thus to provide the 
symmetric JTD spin-phonon coupling. In Sec. 3 we explore the spin-phonon interaction 
assuming a thermodynamical limit. By using the well known technique of the Holstein- 
Primakoff representation we derive an analytical result for the amount of distortion 
and the spin ordering of the system. The experimental requirements for the physical 
implementation and readout of the final state of our model are discussed in Sec. 4. 
Finally, in Sec. 5 we conclude our findings and present further interesting phenomena 
to be explored. 
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2. The Trapped-Ion EgDe Jahn- Teller-Dicke model 
2.1. Trapped-ion radial vibrational Hamiltonian 

We consider a crystal of identical ions with mass M and charge q confined in a linear 
Paul trap along the z axis. Each ion has two metastable internal levels with energy 
separation uq. The system is described by the Hamiltonian (/5 = x,y,z and h = 1 from 
now on) 

TV „ N 2 

i=l /3 i=l 

The first term in Hq describes the energy of the two-level systems with af being the 
Pauli matrices for ion i. H^i\^ is the vibrational Hamiltonian, which contains the ions' 
kinetic energy and the potential energy of the ion crystal. The latter consists of the 
effective harmonic potential and the mutual Coulomb repulsion 

P i=l i>j ' * ■'I 

where is the position vector operator of ion i and Uf^ denote the trapping frequencies. 
In this work we consider the symmetry condition which can be achieved by 

proper adjustment of the trapping voltages or by special design of the trap geometry [34]. 
For sufficiently strong radial confinement (aJr ^ Uz) ions occupied equilibrium positions 
^io = (0, 0, z^) along the z axis. The latter are determined by the balance between the 
Coulomb repulsion and the harmonic trapping force, which yields {dV/dfi)f^-=ff-^ = 0. 
At low temperature the ions undergo only small oscillations around the equilibrium 
positions, namely 

ri = Sr^^iCx + dry^iCy + (z° + 6rz,i)ez, (3) 

where (5r/3^j are the displacement operators. 

The radial vibrational spectrum is essential for the implementation of our idea. 
For that reason we discuss here its main characteristics, a more complete discussion 
can be found, for example in [35,36]. First we notice that a suitable length scale 
is given by d^ = / {Mul). Accordingly, we define a dimensionless equilibrium 
positions, 2;'° = z^/d^. Making a Taylor expansion of the potential (2) around 2;'° 
and neglecting (Jr^ j(5r^ j, 5r^ ^ and higher order terms, the radial vibration is described 
by the Hamiltonian [a = x,y from now on) 

= XI ^ + Yl Yl (4) 

a i=l a i,j=l 

Note that within harmonic approximation of the potential (2) the radial motion is 
decoupled from the axial motion. The collective vibrational frequencies Un = uJt-^/k^ 
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Figure 2. Eigenfrequencies a;„ of the radial collective vibrational spectrum for a linear 
ion crystal with N = 20 ions and w^/wr = 0.1. The highest vibrational frequency 
is the cm mode at uj,-. The frequency splitting between ojc.m and the second-to- 
highest rocking mode Wioc is independent of the number of ions, b) The lowest energy 
vibrational frequency w at as a function of the number of ions for different aspect ratio 
w^/wi- = 0.2, 0.15 and 0.1. As the number of ions N is increased, a critical value occurs 
for which the frequency ujn becomes imaginary, consequently the ion crystal undergoes 
a structural transition into a zigzag phase. The critical N varies for different aspect 
ratio ujz/uJr- 



can found by solving the eigenvalue problem, 

TV 

^/C,,6^^„ = K„6-„, (5) 
i=i 

with the matrix /C given by 

^^J= ( 1 - ^ I ,.0 _\.0 13 ) + ^ U/0 1 ./0|3 - '^-^•) • 

The equilibrium positions for the ions in natural units depend on the number of ions 
only, and thus, the radial vibrational modes are governed solely by the ratio ooz/u)^. 
The vibrational Hamiltonian i^rad can be diagonalized by defining 

N N 

Sra,i = &"„g°(al,„ + Pa,i = i ^ &°nPn(«i,n " (7) 

n=l n=l 

with g° = 1/ y/2Mun, and p° = -\/ Mun/2. By substituting (7) in ifrad we obtain a set 
of uncoupled collective vibrational modes, namely 

^ / . i\ 

(8) 



a n=l ^ 



One important feature of the radial vibrational modes is the fact that the highest 
two eigenvalues do not depend of A^. In particular, if we order the vibrational 
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eigenfrequencies by their energy, we find that the highest energy cm mode ooi = Uc.m = 
ujj- and second-to-highest energy rocking mode 102 = cOmc = y^^? — independent 
of A^, Fig. 2a. This is an important feature, which allows to resolve sidebands to 
the highest frequency Ucm with a resolution independent of A^. Increasing the ion 
crystal undergoes a transition to a zigzag phase, which limits the number of ions one can 
achieve for a given aspect ratio w^/wr [37,38]. In Fig. 2b, we plot the evolution of the 
lowest energy radial mode un with A^, showing the maximum A^ allowed for the linear 
configuration to be stable. In the following we will show that the strong spin-phonon 
coupling in the JTD model induces a structural transition accompanied with magnetic 
ordering, which is a quantum version of the classical zigzag phase. 

2.2. Jahn- Teller E^e spin-phonon coupling 

We describe the interaction of the trapped ion spin ensemble with the collective 
vibrational modes. Those interactions can be induced either by laser dipole forces or 
by magnetic field gradients. Here we focus on the latter technique since it naturally 
implements symmetric couplings in the x-y plane, and it also avoids undesired effects 
like high-order terms in the Lamb-Dicke expansion and the spontaneous emission 
decoherence [28,29]. Let us assume that the ion crystal interacts with an oscillating 
magnetic quadrupole of the form 

Bit; X, y) = bf{t){e^x - Cyy). (9) 

Such a field can be created in a micro-structured planar ion trap, recently experimentally 
demonstrated [30], which contains two wires parallel to the linear ion crystal (Fig. 1) 
[39]. The magnetic field affects only the radial motion of the ion crystal and thus 
the motion along the z axis can be safety neglected. We consider a time modulation 
fit) = (cosz/fet + cosi^rt) to control the couplings. The magnetic dipole interaction is 
described by the interaction Hamiltonian 

N 
i=l 

where /ij = fixO'i + f^yO^^ is the magnetic dipole moment operator of the ion i, and we 
assume the condition = Hy = jji. To control the spin-phonon couplings we choose 
driving frequencies 

t'b.r = (Wo - Wo) ± (Wc.m - t^)- (11) 

The goal is to drive spin-flip transitions with detuning ujq as well as blue- and red- 
sideband transitions of the cm mode cuc.m with detuning iw. The applied bichromatic 
magnetic field in x-y plane, establishes Jaynes-Cummings and anti-Jaynes-Cummings 
interactions, which couple the internal and the motional states of the ions [30,40]. The 
Hamiltonian in the interaction picture with respect to Hq is given by 

h ^ 

Hi= -yY. ^r^.m^t^""''' + ^re-'^"*)/(t) 

i=l 
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h ^ 

- iy E ^rym^t^'""'' - ^re^'^"*)/W, (12) 

i=l 

where af are the Pauh spin-flip operators. The displacement operators 5ra,i are recast 
in terms of collective operators by means of Eq. (7) such that we can in a controlled 
way choose the driving frequencies to pick the radial cm mode as the only resonant 
one. For this, the following set of conditions has to be satisfied, 

A,a;,a;o < t^c.m - Wroc, (13) 

where A = —ixq^hj^^ is the spin-phonon coupling with go = Q'l being the size of the 
cm wave packet. The latter conditions ensures the approximation that any vibrational 
mode but the cm one can be neglected in a rotating wave approximation. Consider as 
an example Zeeman ^^Ca"*" qubits with transition frequency a)o = 20 MHz confined in a 
planar trap with radial frequency = 2 MHz, the first condition in Eq. (13) is justified. 
Assuming Uz/uJv = 0.2 the frequency splitting is approximately (wc.m — i^roc) ~ 40 
kHz. With current ion-trap technology a spin-phonon coupling of order of A ~ 5 kHz 
is achieved by magnetic field gradient 6 = 25 Tm~^, which allowed to neglect the 
contribution of the off- resonant terms in Eq. (12). Under those assumptions we can 
approximate the interaction Hamiltonian by 

Hi = -^(4e'"* + a,e^'"*)(J+e'"o* + e-'"°*J„) 

+ i-^(aje'"* + aye-'"*)(J+e'"»* - J^e"'""*), (14) 

Here and aj^ correspond to the annihilation and creation operators of the cm 
phonon, respectively. Note that the factor N~^^'^ in (14) appears due to the excitation 
of the radial cm modes, wherein the spin-phonon coupling scales as b'ji = N~^^'^. 
Since the ions are equally coupled with the phonons we have introduced the collective 
spin operators J+ = cr^^ (J^ = J_) and = 1/2 ^^^^uf, which describe the 

combined ionic pseudospin of length j = N/2. The collective spin basis is spanned 
by the Dicke states \j,m), which are eigenvectors of J^|j, m) = j{j + l)\j,m) and 
Jz\j,m) = m\j,m), respectively. The Hilbert space of the total system is spanned by 
the states m) ® \nx,ny)}, where \nx,y) is the Fock state with n^^y phonons. After 
performing the time-dependent unitary transformation F = Q^^'ti^^+^y)+^'^otJz ^ such that 
-f^jTD = F'^ HiF — ihF'^dtF, we express the Hamiltonian (14) as 

HjTD = uj{n^+ny)+ujoJz^ — ^{J++J_){al+a^)+i^={J+-J-){al+ay).{15) 

V4j V4j 

Hence we arrive at the realization of the collective JTD model, which describes a two- 
degenerate vibrational modes coupled to the effective spin ensemble by the symmetric 
JT coupling. The Hamiltonian (15) is a multi-particle extension of the E^e model in 
molecular and solid-state physics. The trapped ion realization of JTD model allows for 
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Figure 3. Coupling pattern of the relevant ionic and vibrational states |j, m)|nr,n/) 
for j=l. Due to the symmetries in spin-phonon coupling and parity of the Hamiltonian 
(15), the quantum number n,. — n/ — to is preserved and the Hilbert space is decomposed 
into subspaces of positive (negative) parity. Here we show the non- vanishing couplings 
(solid lines) between states with rir — ni — m = 1 and positive parity. 



easy tuning of the effective spin and phonon frequencies by adjusting the detuning and 
the spin-phonon couphng via the magnetic gradient. 

It is convenient to rewrite the Hamiltonian (15) in terms of right and left chiral 
operators [24] 

4 = ^(4 + i4). ai = ^{ax + iay), (16) 

which can be used to express the z component of the total angular momentum 
Lz = XljLi Lj=nr— ni. Using (16), the Hamiltonian (15) is expressed in the form 

ifjTD = w(4«r + a|a;) + uqJz + ^= J+(4 + ai) + ^= J_(ar + a]), (17) 

which shows that in the JTD model the creation of collective atomic excitation is 
accompanied by the creation (annihilation) of right (left) quantum of angular momentum 
and vice versa. 



2.3. Symmetries 

Due to the symmetry in the spin-phonon interaction the JTD Hamiltonian (15) is 
invariant under the combined application of a rotation in the x-y plane 



ax 




COS 


— sin0 






Ojy 




sin 


cos 




a'y 



and a phase shift J+ — j- e~"^J+. Hence, the JTD model is U(l) invariant, with the charge 
C = Lz — Jz being the group generator, [H, C] = 0. This implies that the Hilbert space 
is decomposed into subspaces with a well-defined quantum number ra,. — ri; — m. Fig. (3). 
Because the Hamiltonian (15) is quadratic in the spin and phonon operators, it is also 
invariant under the application of the parity operator H = exp[i7r + ni + Jz+j)]- The 
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Hilbert space of the total system is thus additionally decomposed into two noninteracting 
subspaces with even and odd number of total excitations [41]. 

2.4- Holstein-Primakoff representation 

In order to study the critical behavior of a collective JTD model (15) in the limit j — )■ oo, 
we use the Holstein-Primakoff transformation, whereby the spin-A^/2 degree of freedom 
is expressed in terms of single mode bosonic operators, namely J+ = — b^b, 

J- = \plj~—Wbb^ and = b^b — j. This transformation preserves the spin algebra and 
allows to convert the JTD Hamiltonian (15) into the Hamiltonian 



which describes three coupled bosonic field modes. This approach is the basis for the 
theoretical discussion in the following section. 

3. Magnetic Structural Phase Transition 

By taking the thermodynamical limit j ^ oo the Hamiltonian (19) can be rewritten as 
follows 



which is bilinear in the bosonic field operators and therefore it can be exactly 
diagonalized by means of a Bogoliubov transformation, which yield (see Appendix A) 




(19) 



ifjTD = uj{nx + Uy) + ujo{b% - j) + —{b^ + 6)(4 + ax) 
+ i^{b^ ~b){a\ + ay), 



(20) 




(21) 



The eigenfrequencies e'^ can be found by solving the eigenvalue problem 



3 




(22) 



for the matrix 
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Figure 4. The energy spectrum of the collective E(X)e JTD model in the limit j oo as 
a function of the spin-phonon coupling A with ^ = 0^0. For A < Ac the eigenfrequencies 
e'p are given as a solution of Eq. (22). At A = Ac the system undergoes a magnetic 
structural phase transition with broken U(l) symmetry. The new eigenfrequencies e'p 
for A > Ac are given by Eq. (28). 



with m± = (1 ±\y^2/ujoUj)^^. The range of vahdity of the eigenfrequencies e'^ is hmited 
for A < a/ ujoUj/2. Indeed, the requirement for hermicity of the matrix (23) is hold for 
A < Ac, with Ac = ^JuJ^p^|2 being the critical coupling. Fig. 4. The vacuum state |0i) 
of Hamiltonian (21) is defined by the condition Cp |0i) = 0. It is straightforward to 
show that the mean-value of the displacement of the cm radial coordinates vanishes, 
(^'"z) = k^'^y) = 0. The phase is characterized with zero phonon excitations and 
collective spin pointing along the z axis, (J^) = — j. 

We may find a simple physical interpretation of the critical spin-phonon coupling 
Ac. Indeed, it is well known that in the presence of spin-orbit coupling, the minima of 
the lower adiabatic potential surface APS (effective nuclei potential in the molecular 
physics) for A < Ac appears at the origin, while for A > Ac APS has a sombrero shape. 
Increasing the spin-phonon coupling, the energy is minimized by breaking some spatial 
symmetry and thus leads to a JT distortion [2]. In order to quantify amount of distortion 
and the spin ordering in the ion crystal above the critical coupling Ac, we follow the 
general procedure introduced by Emary and Brandes in [31] for the quantum Dicke 
model. We displace each of the bosonic modes aj. — )■ aj. + ^/al, a^^ ^ a^y + y^o^, and 
6^ —J- 6^ — a/t*, where ax, o.y, and 7 are generally complex parameters in the order of 
j. By using the Holstein-Primakoff representation and by substituting the displaced 
operators, the Hamiltonian (19) becomes 

ifjTD = uj{a\ax + ^Joi^^x + -Ja'xa^x + |aa;|) + ^(0^^ 

\ay\) ^ u^{\)^h - 
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Figure 5. A magnetic structural phase transition. For spin-phonon coupling A < Ac 
the ions' equilibrium positions are aligned along the z axis, i.e. a;o = yo = and the 
only non-zero projection of the total collective spin operator is ( Jj) ~ —j. Increasing 
A > Ac the system undergoes a magnetic structural phase transition, wherein the radial 
ions' equilibrium positions are displaced by xq = qoy/oiZ and yo = Qo^/c^, accompanied 
by a ferromagnetic spin ordering, {J"^ y) ^ 0. 



X k 



where 



fc = 2j-|7l; 



k 



(24) 



(25) 



The parameters a^^y and 7 can be found from the condition that all terms linear in the 
bosonic field operators in Eq. (24) are canceled 



s2)cos0, = — s^) sin< 



|7| = Vj(1-5), (26) 
with s = A^/A^ and (p = arg(y/7), the second being arbitrary. The latter reflects the 
arbitrariness in the choice of a direction in spontaneous symmetry breaking, Fig. 5. 
Again, in the limit j — > 00, the Hamiltonian (24) can be brought to the diagonal form 
(see Appendix A) 



H 



JTD 



E 

p=2 



1 



u 4A2 



A^ 

2uj 



(l-.).(27) 



The new frequencies solution of the eigenvalue problem 



i=l 



for the matrix 



B"- 



(28) 



(29) 
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Coupling X (units of m) 

Figure 6. The mcan-phonon number {{ux) + {ny))/j as a function of tlie spin-plionon 
coupling A. Tlie numerical results for 10 (red circles) and 20 (blue triangles) ions are 
plotted together with the mean- field solution (black line). For A < Ac the phase is 
characterized with zero mean-phonon number. A position reordering is observed for 
A > Ac accompanied with non-zero mean-phonon number. The magnetic structural 
transition becomes sharper with increasing N. 

with = + l^j, u = (^^^ - ^) and M± = (1 ± s), respectively. In contrast to 
(22), now the frequencies e'^ remain positively defined in the region A > Ac, Fig. 4. 

The mean- value phonon number with respect to the new vacuum state IO2) of the 
Hamiltonian (27) with Vp IO2) = is = ctx,y, indicating a non-zero radial phonon 

excitations. Fig. 6. The collective displacement of the cm mode implies a position 
reordering of the ions' equilibrium positions in the radial x-y plane. Indeed, from Eq. 
(7) it follows that the new radial equilibrium positions are Xq = qo^/a^ and yo = qo^^ya^. 
The structural transition also is accompanied with the ferromagnetic spin ordering, 
( Jj) = j^(l — s^) cos^ (f), (Jy) = j^(l — s^) sin^ (f), respectively. Fig. 7. We note that, the 
magnetic structural transition breaks the continuous U(l) symmetry of the JTD model, 
which refiects to the energy spectrum, namely one of the eigenfrequencies corresponds 
to the gapless Goldstone mode, e" = 0, see Fig. 4. 

4. Preparation and Detection of the Magnetic Structural Phase Transition 

In the following we discuss the implementation of our model in a realistic trapped 
ion experiment. Consider an ion crystal which consists of ^^Ca"*" ions with qubit 
states encoded at the Zeeman 5*1/2 levels, |t) and H). The experimental sequence is 
started by ground state laser cooling of the radial cm modes and pumping spins to 
= \ii ■ ■ ■ D- Because the proposed method for implementation of our model is 
based on magnetic field interaction, the spectator modes could be only Doppler cooled. 
This is a key advantage compared to the laser-ion interaction, where the spin-phonon 
coupling would depend on the spectator modes by the Debye- Waller factor, which is 
a significant source of decoherence [42]. After the preparation of the initial state. 
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Coupling A, (units of co) 

Figure 7. The mean-field result for ((J^) + {Jy))/p as a function of the spin-phonon 
coupling A. The numerical results for 30 (red circles) and 40 (blue triangles) ions are 
compared with the mean-field solution. A creation of macroscopic spin-coherence is 
observed for A > Ac, which is an analog to the super-radiance phase in the Dicke model. 

the coupling A is slowly increased, relatively to the energy gap. For a moderate ion 
crystal with = 10 the energy gap around the critical coupling is approximately 
AE' ^ 0.25ciJ. The adiabatic condition requires T ^ AE'"^, where T is the interaction 
time. Assuming u = ujq = A kHz, the adiabaticity is hold for T ^ ms which is 
comparable with the heating time [30]. The magnetic structural phase transition can 
be detected by measuring either the radial displacement or the spin population at the 
end of the quantum simulation. Assuming spin-phonon coupling A ~ 5 kHz at the end 
of the quantum simulation and direction of spontaneous symmetry breaking </> = 7r/4, 
the equilibrium positions are displaced in the radial x-y directions by xq = ~ 2go- 
Such a structural transition can be detected by laser induced fluorescence, which is 
imaged on a CCD camera. The detection of the magnetic ordering can be performed 
by measuring the spin populations. For particular ion species, an illumination of the 
crystal with resonant light near 397 nm and 866 nm would provide a spin- dependent 
laser fluorescence, namely all spins up emit light and appear bright while spins down 
remain dark. 

5. Conclusion 

We have presented a proposal for the physical realization of the collective JTD model 
based on a linear ion crystal. We have shown that the JTD model exhibits a 
magnetic structural phase transition in the thermodynamical limit. Beyond the critical 
coupling the continuous U(l) symmetry is spontaneously broken which leads to collective 
motional displacement of the radial coordinates and creation of macroscopic spin- 
coherence. The features of the magnetic structural transition can be easily measured 
in the mesoscopic ion crystal by laser induced fluorescence. All parameters can be 
tuned by changing the detuning and the magnetic field gradient. In future we will 
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investigate the JT effects in 2D ion crystals, which are relevant to orbital physics in 
solids. Furthermore, the ion crystal also can serve as a platform for studying non- 
equilibrium phenomena and effects of decoherence in such many-body systems, which 
are computationally intractable. 



Acknowledgments 

This work has been supported by the Bulgarian NSF grants D002-90/08, DMU-03/103, 
DMU-03/107, the EU 7th Framework Programme collaborative project iQIT and RyC 
Contract No. Y200200074. P. A. Ivanov acknowledges the COST Action MP IOTA 
1001. 



Appendix A. Diagonalization of the Mean-Field Hamiltonian 

Appendix A. 1. Normal Phase 

We start with the diagonalization of the Hamiltonian (20). It is convenient to work in the 
position-momentum representation by introducing position and momentum operators 
for each of the bosonic modes, 

^ 1 



X + 



y 




2 ^/2uJ^ 



(Al) 



where the quantum oscillators have frequency u in the x-y plane and Uq in the z 
direction. The transformation gives 



H 



JTD 



[X 



y 



^2 



1 

2 



z' + -{pI + pI + pI) 



+ \\/2uoU}zx — A 



PzPy - eo, 



(A.2) 



with Co = wo(j + 1/2) + ^- The effective Hamiltonian (A.2) describes a system of 
three quantum harmonic oscillators which are coupled through position and momentum 
dependent couplings. To express (A.2) as a set of uncoupled oscillators we need first to 
eliminate the momentum dependent interaction term. To achieve that first we rotate 
the coordinate system along the x axis with the matrix 

1 

(A.3) 







1 1 
^ 

V2 V2 J 



such that {x,y,zY = i?^(gi, g2, gs)^ and {px,Py,PzV = Rx{pi,P2,PzY ■ The transformed 
Hamiltonian is given by 



H 



JTD 



/ 2 2 I 2^ 



1 



^2^2 I Pj ■ 

2^' + 



/ 2 2 

(e ^3 



P3_ 



+ X^/uJ^qi{q3 - ga) 



e )q2q3 - eo, 



(A.4) 
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with = (cj^ + uJq)/2. Hence, the momentum dependent couphng vanishes but as a 
consequence two of the effective harmonic oscillators acquire different effective masses 

m+=(l + \\[^\ , m^=(l-\.[^] . (A.5) 

V V J V V (^o^J 

The Hamiltonian (A. 4) can be rewritten in a compact form as follows 



rr PI 
-njTD — — 



2m-i 



pI 

2m_ 



1 ^ 

9 



2 ^ -^ij (lilj - (^0- 



(A.6) 



Here B^:^^ is (3 x 3) real and symmetric matrix, given by 



(1) 



B^' 



2 _ Jl 



uj — e 

^2 



(A.7) 



—Xy/uiQUj e 

Still, Hamiltonian (A.6) is not in the desirable normal mode form, because the quantum 
oscillators have different effective masses. To overcome this problem we normalize the 
position operators q'l = qi, q'2 = 1/^+^2, q'3 = ^/rnZq^ and momentum operators p'^^ = pi, 
p'2 = P2I p'2, = Ps/ V"^-' respectively, to obtain 

^JTD = 9 J^P? + 9 XI ^'iM - ^0, (A.8) 



i=l 



where 



B' 



m-i- 2m-|- 2^m+m- 



(A.9) 



2-y/m+m._ 2m- 

To find the collective spin-phonon modes, we solve the eigenvalue problem 



(A.IO) 



for the eigenfrequencies e'^ and the eigenvectors v'^^^ with p = 1,2,3. Finally, we 
introduce a new set of bosonic field operators by the relation 

3 



3 /(p) 



q 



1 V^P 



and arrive at the diagonal Hamiltonian 

3 



H 



JTD 



p=i 



clcp + - I - Wo I J + 



(A.ll) 



(A.12) 
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Appendix A. 2. Magnetic Structural Phase Transition 

In order to diagonalize Hamiltonian (24) in the limit j — t- oo we perform the following 
two steps, i) Expand the Hamiltonian (24) as a power series in 1/k and neglect the 
terms in order of j in the denominator, ii) Eliminate the terms in (24) which are linear 
in the bosonic operators by the condition, Eq. (26). The resulting Hamiltonian becomes 



ifjTD = (^{nx + riy) + 7^(1 + s)h^h + -Vl + s x 



(4 + ax){h^ + h)+ i^Vl + s(aj + ay){h'^ - h) 



2^lT~s 

+ sm0(a; + ay)] + — ^^^^ x 

(gi^ftt + e-'*6)2 - go, (A.13) 



with 



eo = j{- + ^} + f(l-s). (A.14) 
oj 4A^ loj 

We can further simplify (A.13) by applying the following transformations 
Ox = dx cos (j) — dy sin 0, 

tty = dx sm (j) + dy cos (f), (A. 15) 

and e^'I'l)^ — )■ 6^. Then, the Hamiltonian reads 

i^jTD = a;(44 + dldy) + ^(1 + s)h'^h 

+ ^7^(&^ + ml + dx) + i^v/r + 7(6t - h){dl + rf,) 
4a; (1 + s) 

Following the standard procedure, we introduce the position and momentum operators 
for each of the bosonic modes 



with Cj = (a;o/2s)(l + s). The Hamiltonian (A. 16) in the position-momentum 
representation is given by 

i/jTD = y (X^ + Y^) + + \{Pl + P^) + ^ 

+ -^PzPx-X\l + s)ZY 
1 + s 

-cu-|-eo, (A.18) 
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with m = (l + s)^/4. It is convenient to normalize the position and momentum operators 
in X and y directions as follows 



\/rnPx Px, 

VrnPy Py, 
Then the Hamiltonian reads 
rl — 



X 



'm 
Y 



Y. 



(A.19) 



~,2 



iJTD 



{X' + Y' 



-^(Pl + + P|) + -PzPx - 2\^mZY 
2m m 



— u 



Co- 



(A.20) 



Similar as before, the diagonalization proceeds by nullifying the momentum dependent 
interaction term in (A.20), which is achieved by rotating the coordinate system along 
the y axis with the matrix 



^ ^ 

1 

^ 

V2 V2 J 



(A.21) 



such that {X,Y,Zf = Ry{Qi,Q2.QzY and {Px,Py,PzV = Ry{Pi. P2. Pzf ■ After 
performing the rotation we obtain 



H 



JTD 



Ql 



PI 



PI 



2mi 2m 2m^ 
V2mX^Q2{Q3 - Ql 



2 m^ 2 
+ + muQiQs 



eo, 



(A.22) 



respectively. Two of the quantum oscillators 



(A.23) 



with e=[i + ^) e.ndu=[e-ii, 
acquire new effective masses 

m m 

with M± = (1 ± s). Again as before to express the Hamiltonian (A. 6) in normal mode 
form we normalize the position and momentum operators, namely: Q'( = ^Jm{Q\^ 
Q2 = \fmQ2 and = ^frniQz and, respectively, P'( = Pxj ^/m{, P'^ = P2I \pm, and 
P3 = P^l ^/m^. Then the Hamiltonian becomes 

3 ^ 3 



eo, 



(A.24) 



4 = 1 



where 



B". 



\^^2M^ 



(A.25) 
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The new eigenfrequencies are obtained by solving the eigenvalue problem, 

E^:>;^^^=4V^^^ (A.26) 

We find that one eigenvalue e" = corresponds to a free mode. The latter is the 
Goldstone mode related to the breaking of the U(l) symmetry. We may define bosonic 
creation and annihilation rp operators for the nonzero energy modes by the 
relation 

p=2 V-^^P p=2 ^ ^ 

Submitting (A.27) in (A. 24) we obtain the following diagonal Hamiltonian, which refers 
to two decoupled oscillators, 

H... ^ ±4 {rW^iy^-y+s)-, g . ^)_^,1_.,.(A.28, 
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